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Abstract 

We perform an analytic calculation of the one-loop amplitude for the VF-boson mediated process 
—7- duQQU retaining the mass for the quark Q. The momentum of each of the massive quarks 
is expressed as the sum of two massless momenta and the corresponding heavy quark spinor is 
expressed as a sum of two massless spinors. Using a special choice for the heavy quark spinors 
we obtain analytic expressions for the one-loop amplitudes which are amenable to fast numerical 
evaluation. The full next-to-leading order (NLO) calculation of hadron -|- hadron — )• VF(— )• ev)hb 
with massive 6-quarks is included in the program MCFM. A comparison is performed with previous 
published work. 
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FIG. 1: Example of a partonic process contributing to Wbh production 
I. INTRODUCTION 

One of the most interesting channels currently under study at the Tevatron and the LHC 
is the final state containing a VT-boson and jets, where some or all of the produced jets are 
tagged as containing a bottom quark. Several interesting partonic processes contribute to 
this final state, for example, 

• W + H{^bb) 

• W + Z{^bb) 

• b + t{^W^ + b) 

• i+t{^W+ + b) 

A background to these processes involving heavy bosons and fermions, is the QCD and 
electroweak process occurring in the collision of hadrons Hi and -ff2, 

Hi + H2^W + b + b + X . (1) 

An example of a partonic subprocess contributing to this process is shown in Fig. 1. Next-to- 
leading corrections to this process were first considered in [1] , working in the approximation 
in which the 6-quark is considered massless. Since in many analyses the 6-quark is required 
to have a minimum pt-, typically 15 GeV or more, in order to be efficiently tagged, the 
neglect of the mass of the bottom quark is expected to be a good approximation. 

After these initial studies, this same process was considered in refs. [2-5] , without making 
the approximation nib = 0. These studies, performed without the inclusion of decay products 



of the W boson effectively confirmed that the nif, — approximation is good to a few percent, 
with the difference parametrically suppressed as m1/p\. Thus retaining the mass of the h- 
quark extends the prediction to lower values of pt- In addition, it allows us treat the case 
where the two 6-quarks end up in the same jet and the case when one 6-quark is either too 
soft or too forward to be tagged. These kinematic configurations can be important for the 
Higgs search at the Tevatron and the LHC [6, 7]. 

The purpose of this paper is to repeat the calculation of ref. [2], i.e. including the effects 
of a finite 6-quark mass, but also including the spin correlations present in the VF-boson 
decay. The calculation is performed using the spinor helicity formalism, with the calculated 
amplitudes represented as analytic formulae. Although the formulae are not compact enough 
to be presented in their entirety here, they do lead to an efficient and numerically stable 
code. 

Following a four-dimensional unitarity based approach [8, 9], we will construct the loga- 
rithmic parts of the virtual amplitude using multiple cuts. Developments to this technique 
utilizing complex momenta [10-13] allow us to compute analytic expressions for the coeffi- 
cients of the known scalar integrals. 

The bulk of our paper is dedicated to a description of the calculation of the one-loop 
corrections to the process duQQii retaining the mass for the quark Q. The detailed 
plan of this paper is as follows. In section II we present our method for dealing with 
massive spinors using spinor helicity techniques. Section III gives a precise definition of 
the amplitude that we wish to calculate, including the decomposition into colour stripped 
amplitudes and the further decomposition into one-loop primitive amplitudes. Sections IV, 
V and VI illustrate our calculation of the leading colour primitive amplitude, A''^, of the 
sub-leading colour primitive, Af, and of the primitive amplitudes containing a closed loop 
of fermions, A^^ and A^^. Section VII presents the renormalization counterterms. After 
describing the implementation of the calculation into MCFM in section VIII and comparing 
with earlier work, we draw some conclusions in section IX. 

II. TREATMENT OF MASSIVE SPINORS 

A method for dealing with a massive particle in the context of the spinor helicity method 
has been given by Kleiss and Stirling [14]. A massive momentum can always be represented 
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as a sum of two massless momenta. Thus the spinor solution for a massive particle can be 
expressed in terms of massless spinors by decomposing the physical momentum in terms 
of the two massless momenta. If we are ultimately going to sum over the spin degrees of 
freedom, the only constraint that the massive spinors must satisfy is that they should give 
the standard result for the spin sum after averaging over polarizations, namely. 



Us{p, m)us{p, m) = ^+ m , 
^Vs{p,m)vs{p,m)^j^-m. (2) 

We can decompose a massive vector into two massless vectors by introducing an arbitrary 
massless reference vector, 77, 

In this equation p^ is a massless vector. The details of our spinor product notation are given 
in Appendix A. The definitions of massive external spinor wave functions are, 

u±{p,m;p\r]) = -^-^^ (77 ^ |(/ + m) , 

v±{p,m;p\rj) = m)|77±) , (4) 

where the subscripts ± label the spin degrees of freedom. In the massless limit these labels 
correspond to the helicity quantum numbers, but in the massive case they have no such 
interpretation. Treating the spinors as independent functions of p^ and 77 i.e. ignoring the 
constraint in Eq. (3) we can show using simple manipulations that, 

[p^ 7]] 



m 



U-{P:m;r],p) = u+{p,m;p ,rj) , 



^^—^v+(p,m;r],p^) — v^(p,m;p^,r]) . (5) 
m 

This has the attractive feature that amplitudes with different spin labels can be obtained 
from one another by exchanging p" and 1]. This method has been used in the calculation of 
one-loop corrections to top production [15, 16]. 



A. Special choice for massive spinors 



In this paper, however, we adopt a different approach. By making a specific choice for 
the vector rj in terms of other vectors in the problem we can simplify the calculation of 
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individual amplitudes. In addition, we shall find that for our particular choice of the vector 
77, one-loop results for the colour suppressed primitive amplitude can be obtained directly 
from the corresponding massless amplitude. 

The implementation of the Kleiss-Stirling scheme appropriate for the case where we have 
pairs of massive particles is due to Rodrigo [17]. The two massive momenta p2 and with 
equal masses, corresponding to the momenta of the antiquark and quark respectively, are 
written in terms of lightlike momenta k2 and k^, 

M _ 1 + /3 , ^ _L 1 - /3 , M 

P^ = ^fc3" + ^fc2^ (6) 



where. 



= Vl - 4mVs23 , ^± = ^(1±^), (7) 
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and S23 = (P2 + Pa)^ = 2A;2 • k^. The decomposition of Eq. (6) has the advantage that 
momentum conservation is preserved, P2+P3 — A;2 + ^3- The inverse transformation is given 

by, 

*? = ^K-^p5. (8) 

In the rest of this paper we shall denote massive vectors by Pi, (pf 7^ 0) and massless vectors 
by ki, {kf = 0). For the massless vectors ki we shall further define massless spinors, 

u-{ki) = = \i] , u+{ki) = \i+) = \i) , 
u-{ki) = {i -\ = {i\ , u+{ki) = + \ = [i\ ■ (9) 

In terms of the two massless momenta in Eq. (8) the explicit results for the massive solutions 

of the Dirac equation are, 

^-1/2 ^-1/2 
u±{p3,m) = -^^^^{2'f\{^3 + m) , v±{p2,m) = ^0^{^2-m)\3^) . (10) 



This corresponds to choosing. 



V2 = ^^3, P2 = P+k2 , 

Vs = h, p\ = ^+k3 , (11) 
5 



in Eqs. (3,4). 

With this choice, the results for the massive quark current with spin labels {+, — } and 
{— ,+} have the same form as they would have in the massless limit (i.e. in the limit in 
which Pi — )■ ki for i = 2, 3), 

5'^(3±,2J) = u±{p^,m)Yv^{p2,m) = u±{h)rv^{k2) = {S^hT^} • (12) 

The results for spin labels {— , — } and {+, +} on the heavy quark current are, 

S''iS-,2^) = u4ps,m)Yv-ip2,m) = 2jV__ {k^ - k^Y , 

>S'^(3+ 2t) = tZ+(p3,m)7V(P2,m) = {k2 - hY , (13) 

where the overall normalization is given by. 

Contracting these equations with a Dirac matrix we obtain, 

5^^(3+2^)® {7,} = 2{|2)[3| + |3](2|} , (15) 

S^{3-,2±)^{j,} = 2{|3)[2| + |2](3|} , (16) 

5'^(3q,2^)®{7,} = 2A/-__{|2)[2|-|3)[3| + |2](2|-|3](3|} , (17) 

5'^(3;,25)^{7,} = 2A/;+{|2)[2|-|3)[3| + |2](2|-|3](3|} . (18) 

As an example of the use of these spinors, in Appendix B 1 we outline the calculation of the 
tree-level amplitudes for gQQg scattering, which will appear later as an ingredient in the 
calculation of our one-loop amplitudes. 



III. SETUP 



We shall consider the process, 

^ q{ki) + Q{p2) + Qips) + q{h) + + i{k&) , (19) 

both at tree level and including the one-loop QCD corrections. The process is mediated 
by the exchange of a VF-boson which decays into an antilepton-lepton pair with momenta 
labelled by k^ and ke, as shown in Fig. 2. The V-A structure of the charged weak interaction 
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FIG. 2: Feynman graphs that enter the calculation of the lowest order amplitude. The massive 
quarks are represented by the heavy (red) line, the wavy line denotes a W boson and as usual the 
helical line denotes a gluon. 



ensures that the lepton and massless-quark lines will have fixed helicity. Thus the outgoing 
lepton (6) will be always left-handed and the outgoing antilepton (5) always right-handed. 
With this understanding we can often drop the lepton labels from the specification of the 
amplitude. 



A. Colour decomposition 

As noted above we shall suppress the labels for the outgoing leptons 5 and 6 which play 
no role in the colour decomposition. The colour decomposition for the tree graphs is, 

ATil,, 2q, 3q, Ag) = Vw{s,e) A'r{l„ 2q, 3q, 4,-) (^/^ 6^^ - j^^^ 6^^ ) . (20) 

In Eq. (20) we have introduced the strong coupling constant, g, and the weak coupling, gw, 
defined through. 

The Breit-Wigner factor Vwis) is given by. 

The indices ji, js, ( j^? J a) denote the colour labels of the corresponding quark (antiquark) 
lines. 

At one loop the colour decomposition is given by, 

^6"'°°P(lg, 2q, 3q, 4q) = g^g^cr Vw{s56) 

5]^ 5]^ A6;i(l„ 2q, 3q, 4,-) + 5]^ d]^ A,,2{1,, 2q, 3q, 4,-)] , (23) 
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where the overall factor cr is, 



1 r(l + e)r2(l-e) 



- (47r)2- r(l - 26) ■ ^^^^ 

The interference of the one-loop amplitude with the lowest order, summed over initial and 
final colours, is given by. 



J2 [-^^AInlo = 2gl^g'cr{N^,-l)N,\Vw{s56)\'' 

colours 

X Re{K^^(l„2,3,4,-)]M6.i(l„2,3,4,-)}. (25) 
Therefore Aq-2 plays no role in the calculation of NLO corrections. 



B. Tree level amplitudes 

The tree level amplitudes can be calculated using the diagrams of Fig. 2. The { — h} 
amplitude, expressed in terms of the momenta ki, is identical to the massless result as a 
consequence of Eq. (12), 

;i3)[4 5](6|(l + 3)|2] [42](16)[5|(2 + 4)|3)- 



-^4'^^ni;,2t,3-,4+,5+,67) = 



523=^56=^123 

(13) [4 5] (6|(1 + 3)|2] | 
S23S56S123 J 



S23S56'S234 

flip 



where we have introduced the symmetry, 

flip: (1^4), (2^3), (5^6), []^( 

The { } amplitude is given by, 

- iAT%l- 25, 3-, 4+ 5+ 67) = 



(26) 



(27) 



X 



^^^^^■ (6|(1 + 2)|3](31)-(6|(1 + 3)|2](21) -^ ^ <( flip 



, (28) 



•§23*56*123 

with the mass-dependent normalization factor given in Eq. (14). The amplitudes for 
^*6'''''(l",2-,3+,4+,6-,5+) and ^^'^'^^(l", 2+, 3+, 4+, 6", 5+) can be easily obtained by ref- 
erence to Eqs. (12,13), 



^T^'i^q ' 2q, 3+, 4t, 5^ , 6^ 



/itrcc/i- 0+ 9- A+ p:+ 
6 y^q^'^Q^'^Q^ 7^ i 



Ar(l,",25,3j,4+,5+,67) = -flip 



Atreez-i- r,- o- a+ r:+ c- 



(29) 
(30) 
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FIG. 3: Parent diagrams for the leading colour primitive Ag(l,2,3, 4) and crossed box primitive 
At {1,2,3, A). 

C. Decomposition of one-loop amplitudes into primitive amplitudes 

Following Ref. [18] we can decompose the full one-loop amplitude into gauge invariant 
primitive amplitudes, 

A;i(1„2q,3q,45) = ^1^(1,2,3,4)- J^«(l,2,3,4) + A|?(l,2,3,4)) 

- ^^Alil, 2, 3, 4) - ^Ajf (1, 2, 3, 4) - ^A^^(l, 2, 3, 4) , 

A6;2(1„2q,3q,4,-) = <(1,2,3,4) + ^(A|?(1,2,3,4) + <(1,2,3,4)) 

+ -^^Afil, 2, 3, 4) + ^<(1, 2, 3, 4) + ^<(1, 2, 3, 4) . (31) 

A^Q is the primitive amplitude containing the leading colour box and Af" is the primitive 
amplitude containing the crossed box. The parent diagrams are shown in Fig. 3. These two 
primitive amplitudes are related by. 



Ach/^hi 9/12 q/13 AhA\ _ /ilc/i/ii 0^2 Ahi\ 



Q ' Q ' 9 



(32) 



is the primitive amplitude containing the subleading colour boxes, see Fig. 4, and the 
primitive amplitude for the fermion-loop diagrams (light and heavy), y4g,y4g^ is shown in 
Fig. 5. We note that, just as at tree level, parity ensures that we need only compute a 
reduced set of amplitudes. For A^^ at one loop it is sufficient to calculate three of the four 
possible spin label combinations since (c.f. Eq. 30), 



4'(l,"'2t 3+4+5+67) = -flip 



^6 (Ig ' 3q, 4t, 5^ , 6^ 



(33) 
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FIG. 4: Parent diagrams for the subleading colour primitive amplitude, Af[l,2,3,4:). 
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FIG. 5: Diagrams for the fermion loop primitives (1,2, 3, 4) and ^g^(l, 2, 3, 4). 
D. Structure of the calculation 

In the presence of propagators with vanishing masses on internal lines the one-loop am- 
plitude contains infrared and collinear divergences. In addition, the amplitude contains 
ultraviolet divergences. We regulate all of these divergences using dimensional regulariza- 
tion, continuing the loop integration to = 4 — 2e dimensions. The divergences then appear 
as poles in e. For the primitive amphtudes the divergence structure is quite simple and we 
can separate the amplitude as follows, 



^6 



(34) 



where contains all the divergent pieces. In this equation Aq is any of the primitive 
amplitudes and the index j runs over {Ic, cb, si. If , hf}. As we shall see later the simplicity 
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FIG. 6: Schematic diagram illustrating the decomposition in scalar integrals. 



of the pole structure gives useful constraints on the different terms contributing to the 
amplitude. 

The calculation proceeds by noting that the primitive amplitudes can be expressed as a 
sum of scalar box, triangle, bubble and tadpole integrals and a rational part. 



where the label j runs over the values {Ic, cb, si, If, hf}. In this equation, illustrated in 
Fig. 6, the indices ii,i2,i3 run over all partitions of the cyclically ordered momenta. The 
coefficients d, c, b, a of the scalar integrals in Eq. (35) also depend on the primitive amplitude 
index j. In the following we shall suppress this dependence, but the primitive amplitude 
to which we are referring will be clear from the context. The functions ^40, Bq, Cq, Dq are 
the scalar integrals defined in Appendix C. The mass labels have been suppressed in the 
scalar integrals in Eq. (35) because they are determined automatically by the propagators 
connecting the two external massive quarks. Values for these scalar integrals are obtained 
using the QCDLoop program [19]. By using four-dimensional unitarity methods [10-13, 20], 
we obtain the coefficients, a, 6, c, d in the four-dimensional helicity (FDH) scheme [21]. The 
rational terms require information from beyond four dimensions. 



-^4(l-,2'^^3'^^4+,5+,6-) - di,\i,\i,{hs, h^) Doipi„Pi„Pi,) 



+ J^^*il^2(^3,^2) C'o(Pn,Pi2) ^6ii(/l3,^2) Bo{pi^) + a{h3, h2) Aq 



+ i?^■(l-,2'*^3'^^4+,5+,6-) , 



(35) 
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Scalar integral 


Coefficient 


Scalar integral 


Coefficient 


1. Do{ki,p2,P3;0,0,m,0) 

2. i:)o(^i,P2,P34;0,0,m,0) 

3. Do{ki,p23,k4;0,0,0,0) 


C^l|2|3 
C^l|2|34 
<^1|23|4 


4. Do{p2,P3,k4;0,m,0,0) 

5. Do{pi2,P3,k4;0,m,0,0) 


C^2|3|4 
<^12|3|4 



TABLE I: Scalar box integrals appearing in the leading colour amplitude Agl,2,3,4) and the 
notation used in the text to denote their coefficients. 

IV. THE RESULTS FOR PRIMITIVE AMPLITUDE Ajf 
A. Divergent parts 

For the leading colour primitive amplitude the divergent term which enters in Eq. (34) is 
given by [22], 




B. Calculation of the box coefficients 

The five scalar box integrals that enter the leading colour amplitude are enumerated in 
Table I. Of these, we directly calculate the coefficients of boxes 1-3 and obtain the remaining 
two by symmetry: 

C^2|3|4(^3, ^2) = -flip [(ii|2|3(-^2, -/i3)] , 
C^12|3|4(^3, ^2) = -flip [c?l|2|34(-^2, -^3)] , (37) 

where /is and /i2 represent the spin labels of the quark and antiquark respectively and the 
operation flip is deflned in Eq. (27). When using the Rodrigo choice to decompose p2 and 
P3, box 3 only gives a non-zero contribution for the spin-labels { — , — } and {+, +}. 

The coefficients are computed using quadruple cuts with complex momenta [10]. For 
boxes 1 and 2 the presence of a massive propagator leads to a more complicated parametriza- 
tion of the loop momentum than in the purely massless case. The particular parameteriza- 
tions that we have used are spelled out below. 
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FIG. 7: Diagram for the calculation of the coefficient of the scalar integral Do(ki,p2,P3; 0, 0, m, 0). 
The massive quark is represented by the heavy (red) line and unitarity cuts are represented by a 
dashed line. 

1. Calculation of di\2\3 

The cuts used to isolate this coefficient are depicted in Fig. 7, where we remind the reader 



that we work in the case that kf = kf = and P2 
to the following constraints, 

/2 = o, {l-k^f = 0, H-k.-p^f- 



p| = m?. The loop momentum is subject 



{l-h-p2- pzf 



0, 



(38) 



Denoting li = I — ki we have /i ■ fci = so we must have either ~ |1) or ~ |1], 
depending on the helicity of the internal gluons, in such a way that the 3-point qqg vertex 
is non-vanishing. 

In the case [/i | ~ [1 1 we can write 

[i|7'^;2f|i] _ [i|7'';2f|i] 



'123 



= - A;f 
I1-P2- P3 



2 im+mi] 2 
2 



(39) 



[ll^^ll] 

where the boldface momenta denote massive vectors. By inspection, it is clear that this 
parametrization of the loop momentum satisfies the constraints in Eq. (38). Correspond- 
ingly in the case ~ (1| we can write 

S23 mrii) _ 



'123 



l^l = F - fcf 
K-P2- Ps 



2 (i|;2f(;2 + ^)|i) 2 
2 



(40) 
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Note that the denominator in Eq. (40) simphfies when expressed in terms of massless mo- 
menta, 

(1|^^|1) = /3(12)[2 3](31) , (41) 

As an example, we will calculate the coefficient of this box integral for the spin label 
choice {—,+}■ The box coefficient is given by the formula. 



^i|2|3(-,+) = lj2dlmt'\-,+) 

ha,hb 

^IYI ^3(1-, -/^/^") X A5(/-,4+,-4^3,5±,6-) 

xA,{-l-'^-,2±3- X i-t) [{l-k,-p2f-m'] , (42) 



2 



with the loop momenta fixed according to the constraints above, i.e. l^^ is given by cither 
Eq. (39) or (40) above. In this formula we have for brevity suppressed the gluon particle 
labels. Note that we have not explicitly cut the heavy quark propagator but, equivalently, 
have simply written the full gQQg amplitude and multiplied by the cut propagator. Since 
the gQQg amplitudes with the same hehcity gluons vanish (c.f. Eq. (B13)), the coefficient 
receives only two contributions. 

Let us first inspect the assignment ha = —hf, = +1. We note that it is simplest to 
manipulate the gQQg amplitude, (Eq. (B13) first line), using the Schoutcn identity, Eq. (A6), 
into an alternative form in order to take full advantage of the vanishing of the cut propagator, 

[2 ^123]^ (^1 1^1 ^123] 



-iA^-li, 2+ 3^, /+3) [{h - P2f - m^] 



[kh2z] [2 3] 

[2/i23]'([^i^i23] {k\m-{h\m [SM) 
[ii3][ZiZi23][2 3] 

[2^123]' fn ,^,01 [3^123] 



((,1,131 - I!M|Jpl) ,43) 



[h 3] [2 3] 

On shell, (^i — ^2)^ — = (/i|^|/i] = and we can simply discard the second term in this 
equation. Using results for the other amplitudes presented in Appendix B we find that, 

^1^1^^ [11] [4il23][ni23][6 5] m ' ^^^^ 

where we have simplified the gQQg amplitude further using the decomposition of p2 in 
terms of k2 and k:^. Since the propagator I does not have a simple decomposition in terms of 
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external momenta, it is simplest to eliminate it by multiplying this expression in numerator 
and denominator by a factor {A I) (lusl). After some simplification this yields, 



^^ _ /^+ [4 5f (41) (/i23l) [l/i] [2/i23]^(2/i) 
^1^1^^ ^123 [6 5] (4 /i) [4/123] [/i 3] ' ^ ^ 



This is now in a form where we can use the loop momenta definitions given in Eq. (40) with 
spinors identified by, for example, using the identity, — We can thus make the 

replacements: 

(/i23 1) (i|;2f^|i) - ' ' 

[4M ^[4|(^+^)|1) = (1|2 + 3|4], [2Zi23]^[2|(^ + ^)|l) = [23](31) , ^^^^ 

{2h) ^(21), (4Zi)^(41), 

This results in the final expression, 

_ (13)^[23][45]^(l|;2f|l] 
^1^1=^^ .123(112 + 3|4] [5 6] • ^^'^ 

Adding in the contribution from the other helicity configuration [ha = —hj, = — 1), computed 
in a similar fashion, yields the full result for this coefficient. 



.23(i|;2f|i] 

Ctl|2|3(-,+j - 



[23] /^(6|;2f|l]^(16)[12]V (13)^ [45]^ 



(5 6)(4|2 + 3|1] V/5[13] [23] J (2 3) [5 6] (1|2 + 3(4] 

(48) 



This is in agreement with the coefficient presented in Ref. [18] in the limit that the heavy 
quark mass is taken to zero. 

2. Calculation o/di|2|34 

This coefficient is obtained from the cuts shown in Fig. 8, where the loop momentum is 
subject to the conditions, 

/2 = 0, {l-hf^O, {l-h-p2f^m\ /?234 = (^-fci-P2-P3-M' = 0. (49) 

An explicit solution for li — I — ki is given by, 

- 2 (i|;2f(^ + ^|i)' ^'"^ 
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FIG. 8: Diagram for the calculation of the coefficient of the scalar integral i5o(^iiP2,P34; 0, 0, m, 0). 
The massive quark is represented by the heavy (red) line and unitarity cuts are represented by a 
dashed line. 

which clearly satisfies the constraints in Eq. (49). From this we can also derive a useful form 
for the momentum /1234, 



'■1234 



h-P2-P3-h = -- 



(51) 



2 (l|;2f(^ + ^|l) 

The other solution, corresponding to ~ [1| can be obtained from Eq. (50) by complex 
conjugation. We note that both of these parameterizations can be obtained from Eqs. (39) 
and (40) by simply making the replacements, ^ — t- (^ + ^ and S23 — )■ S234. From these forms 
it is straightforward to compute the coefficients of this box. 



3. Calculation of (1^2314: 

The cuts used to isolate this coefficient are depicted in Fig. 9. This box with massless 
internal lines and two opposite massive external lines is often referred to as the easy box [23] . 
The heavy quark does not participate in the loop, and the loop momentum is subject to the 
following constraints, 

/2 = o, {i-kif = o, (i-ki-p2-P3f = 0, {i-h-p2-P3-hf = 0, (52) 

The explicit solutions for li = I — ki and /123 = i — ki — p2 — ps are given either by, 

l(4|(^ + ^)7ll) 



2 (41) 



_ 1 (417^^+^)11) 

^123 - 2 (41) ' ^ ^ 
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FIG. 9: Diagram for the calculation of the coefficient of the scalar integral Do{ki,p23, k^; 0, 0, 0, 0). 
The massive quark is represented by the heavy (red) line and unitarity cuts are represented by a 
dashed line. 



or by, 



n 



'■123 



2 [41] 
2 [41] 



(54) 



Both of these clearly satisfy the constraints in Eq. (52). 

As an example, we will calculate the coefficient of this box integral for the spin label 
choice { — ,+}. For this spin label choice it turns out that the box coefficient vanishes. The 
argument is quite general and may be useful in explaining the absence of the easy box in 
other contexts, e.g. ref. [24]. The result for the box coefficient is calculated from. 



ha,hb 



X ^3 (^1234 (7 5 4^, ^123) X li "5 2^, 3„, /]_23'' 



(55) 



with the loop momenta fixed according to the constraints above, i.e. is given by either 
Eq. (53) or (54) above. From Eq. (B13) we see that A4(— /j^'^", 2^, 3~, /^23'') vanishes for 
ha = hb- Therefore we need only consider the case ha = —hb. For definiteness let us consider 
the case ha = 1, hb = —1- For this case we find that. 



^3(1-, -/%,/+) 



[Ih 



^3(^1234 g' ^123) 



[1 /] [h 1] 

(^1234 ^123) 
(4/123) (4/1234) 



1/123] ~ |4] . 



(56) 
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Scalar integral 


Coefficient 


Scalar integral 


Coefficient 


1. 


Cn(p2s, A:4;0,0, 0) 


C23|4 


6. Co(A;i,p23, 0,0,0) 


Cl|23 


2. 


Co(Pi2,P3;0,m,0) 


^12 1 3 


7- C'o(j'3)^4)"T')0, 0) 




3. 


Co(pi2,P34;0, m, 0) 


Cl2|34 


8. Co(A;i,p2,0,0,m) 


Cl|2 


4. 


Co(P2,P3;0, m, 0) 


C2|3 


9. Co(A;i,p234, 0,0,0) 


Cl|234 


5. 


Co(P2,P34;0,m,0) 


C2|34 


10. Co(pi23,fe4, 0,0,0) 


Cl23|4 



TABLE II: Scalar triangle integrals appearing in the leading colour amplitude A|f(l,2,3,4) and 
the notation used in the text to denote their coefficients. 

However we see that the conditions ~ |1), I/123] ~ |4] in Eq. (56) are not compatible with 
the kinematic constraints, Eqs. (53,54) which require either, 

~ |1) and 1/123) ~ |4) , 
or |/i] ~ |1] and I/123] ~ |4] , (57) 

respectively. The argument follows identically for the contribution ha = —l,hh = +1. 
Therefore for this external spin label choice {— , +} the box coefficient is zero. For the spin 
label choice {— , — } there is a non-zero easy box coefficient because 744(— /]f'*", 2^, 3^, /^gs*") 
no longer vanishes for ha — hb. 

C. Calculation of triangle coefficients 

The leading colour amplitude receives contributions from the ten triangle scalar integrals 
listed in Table II. The coefficients of Triangles 1-4 are calculated directly as detailed in the 
subsection below. Triangle 5 can be simply obtained by symmetry, 

C2|34(/i3,/i2) = -flip [Ci2|3(-/i2,-/i3)] . (58) 

The coefficients of the five remaining triangles are then uniquely determined by the known 
divergence structure [22] of the amplitude, Eq. (36). The terms in Eq. (36) proportional to 
log(— 2A;i.p2)/e and Iog(— 2A;4.p3)/e fix two of the triangle coefficients, 

c,2 = -f(i|;?|i](-ar) + ^ + ^) , 

V -523 S234 / 

C3, = - ((41^141 H^r) + ^ + ^) . (59) 

V ^23 Sl23 / 
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FIG. 10: Triangle diagram showing the momentum parameterization. All momenta are outgoing, 
Ki + K2 + K3 = 0. 

while the absence of single poles of the form log(s)/e for s e {523,51234,5123,5234} requires 
the following relations, 



Ci|23 = (1|2 + 3|1] 



'^1|2|J 



+ 



d 



2|3|4 



C23|4 



2c?l|23|4 



523(l|;2f|l] 523(4|^|4] (4|2 + 3|4] (1|2 + 3|4](4|2 + 3|1] 

Cl|23 



Ci|234 = (l|2 + 3 + 4|l](^^+ ^^1^1- 



ci23|4 = (4|l + 2 + 3|4] 



523(1|;?|1] 5234(l|;2f|l] (l|2 + 3|l] 

^^121:5,1 , <'^l,2,:-!4 '"l,2:-!l 



(60) 



+ 



2f/i, 2:5,1 



.5i23(4|^|4] 5234(1|?'|1] (l|2 + 3 + 4|l] (1 12 + 3|4] (4|2 + 3| 1] 

These expressions are written in terms of already-calculated box and triangle coefficients, 
the leading order amplitude, and the other triangle coefficients discussed below. 



1. Forde method for triangle coefficients 

We will calculate the coefficients of the triangle integrals using the method of Forde [13] . 
Triangles 2, 3 and 4 will require a slight extension of the formahsm to include one of the 
internal propagators with a mass. We first review the case of three massless internal mo- 
menta, shown in Fig. 10 in order to introduce our notation which differs from that of Forde. 
Defining li and I2 as follows, 

/r = /o^-i^r, it^ = i^,+K^, (61) 

the cut loop momenta (/^ = 0), i = 0, 1, 2 may be written in the following general form, 

C = x,K',^^ + y,KY + \{K\mK'^ + '^{K\\r\K'^ ■ (62) 
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All momenta can be expanded in terms of massless momenta, and K2, 

7 

K2 ^ K^2 + —K\ , 

7 

Ks = -(l + ^)K{-(l + ^)Kl (63) 

7 7 

where 5'^ = Kf and 7 = = • Xj. The inverse relations are, 

_ - (^i/7)i^2" _ - {S2h)K^ . . 

The dot product of Ki and ^2 in Eq. (63) produces a quadratic equation for 7, the solutions 
of which express 7 in terms of the external momenta, 

7± = (Ki ■K2)±VA, a = (Xi • - SiS2 . (65) 

Prom Eq. (62) the massless vectors li can be expressed as a linear combination of the spinor 
solutions for the vectors K\ and K2, 

{k\ = t{Kl\ + y,{K'2\ , [h\ = ^[Kl\ + [K',\ . (66) 

The on-shell conditions = for i = 0, 1, 2 allow us to derive the coefficients, Xi and y^, 

^ M7+^ ^_S2b±Sil 

^1 -Si>S2(7 + '5i) 7(7 + -^2) 

7 7(7^ - S1S2) T - S1S2 

7(7 + ^1) ^^2 S^S2 (7 + ^2) 
1/2 = 2/0 + 1 = ^ X2 = Xo + — = -— . 67) 

Y - S1S2 7 7(7^ - -S'15'2) 

The spinor products can be expressed as follows, 

[^^i] = ^-^[K'2K\] = -][K'2Kl], 



(Ik) 



{Ik) 
[hk] 



Xi - 


Xo 


t 




t{yi- 


-2/1 


X2 - 


Xq 


t 




t{y2 - 


-2/1 


X2 - 


Xi 



7 



2] = = ^[i^^i^^], 



{h h) = ^(2/2 - 2/1) (i^t i^^> = ^ (1 + 7) (^1 ■ 



(68) 
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FIG. 11: Diagrams for the calculation of the coefficients of the scalar integrals Co(pi2,P3; 0, m, 0), 
Co{pi2,P34:',0,m,0) and Co(p2iP3; 0, m, 0). The massive quark is represented by the heavy (red) 
line. The momenta on the external lines are all outgoing. 



Turning now to the case where the propagator with momentum I2 has a mass, we have 



that Iq = 11 = and /| = m^. Since I2 is no longer massless it does not have an expansion 



of the form Eq. (66). The results for the coefficients in the expansion defined in Eq. (66) for 
Iq and 1 1 are, 

Si {7 + S2- vr?) jm^ 



1/1 = 1/0- 



(7' - ^1^2) 
Si 5i (52 (1 + (5i/7)) - m2) 



1 



72 - SIS2 



Xi = Xq — 1 



S2 (7 + Si) 

(7' - ^1^2) ' 

- 7(5*2 + 7) 



72 - S1S2 



■ (69) 



With these results in hand we can then compute the coefficients for triangles 2-4 with 
the loop momentum assignments as shown in Figure 11 by making the replacements: 

• Triangle 2: Ki = -pus, K2 = P3. 

• Triangle 3: Ki = -J91234, K2 = pu- 

• Triangle 4: Ki = -P23, K2 = ps- 

For the case of triangles 2 and 4, we see that the coefficients are particularly simple since 

52 = m\ 

We shall sketch the calculation of the coefficient of triangle 4 because it is the simplest. 
We have 5i = P23 = -523 and S2 = pi = m? and hence see that, 

7± = i^i ■ i^2 ± v^(^r^^F^5^ = -P^S23 . (70) 
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Scalar integral 


Coefficient 


Scalar integral 


Coefficient 


1. 


5o(pi234 : 0,0) 


^1234 


4. 


Bo{pi2;0,m) 


bi2 


2. 


S0(P23;0,0) 


^23 


5. 




b34 


3. 




hr2,, 


6. 




h2:n 


Bo{p2;0,m) 


b2 


Ao{m) 


a 



TABLE III: Scalar bubble and tadpole integrals appearing in the leading colour amplitude 
(1,2,3,4) and the notation used in the text to denote their coefficients. 

The lightlike momenta, Eq. (64), for the two solutions for 7 reduce to. 



Kl^-k2, Kl 



ks , for 7 = 7_ = -P+S23, yo = xo = 



-S23, yo = a;o = +^ ■ (71) 



Kl^-ks, K'^^/3_k2, for 7 = 7+ = 

Using these assignments we find that the formulae for the coefficients of this triangle are 
relatively compact. For example, for the {— , +} spin labels the result is, 

1 



C2|3(-,+J 



2^56 



[45]/ /3_ (13)^(61(1 + 2)13] _ /3^(12) [1 2]^ (6|(1 + 2)|3] 
S123 V (1 2) p [1 3f 



(1 2) [1 2] (6|(1 + 3)|2] (1 3) [1 2] (6|(1 + 2)|3] 



/3[13] 



(13)(6|(1 + 3)|2]) 



/3[13] 

/3_ (13) (16) [4 5] 



+ 



(12) 



/32 (23)(4|(2 + 3)|5](6|(1 + 2)|3] , /3+(3[(2 + 4)|5](6|(1 + 2)|3] 



+ 



/3 (2 4)2 [13] ' (2 4) [13] 

D. Calculation of bubble and tadpole coefficients 



- <^ flip 



• (72) 



The bubble integrals present in the leading colour amplitude are shown in Table III. The 
coefficients of bubbles 1-4 are computed by using the method of spinor integration [11, 12, 
20]. This is straightforward to apply for bubbles 1-3 but requires a small modification for 
bubble 4 due to the effect of the massive propagator. The modified method is described 
in the subsection below. Coefficients of the bubble integrals 5 and 6 are then obtained by 
symmetry. 



^34(^3, ^2) = -fiip[6i2(-/i2,-^3)] , 
^234(^3, ^2) = -fiip[6i23(-^2,-^3)] , 



(73) 
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All occurrences of So(p3;0,m) have been replaced by Bo{p2;0,m.), so we only need to 
determine the coefficient 62- A linear combination of the coefficient 62 and a is determined 
from the known form of the single pole in e, (i.e. the single pole with no associated logarithm 
in the expansion of Eq. (36)), 

g 

62 + m^a = -{-i)Al'''^ - b^s - ^234 - &23 - ^1234 - ^12 - &34 ■ (74) 
Eq. (74) is sufficient to fix the poles and the logarithms, but since, 

Bo(P2;0,m) = (^)'(i + 2), (76) 

it leaves the constant term undetermined. In order to separately fix the coefficients a and 
62 we perform a direct Feynman diagram computation of the tadpole coefficient a to find, 

«(-,+) = 0, 

a(+, -) = , 
mM-,-) = ^(-i)Ar(l,-,2^,3-,4+,5+,67), 
mM+,+) = ^(-^)^r(V,2t,3+ 4+,5+,67) . (76) 

These results for a are valid in any covariant gauge. 



1. Bubble integral with one massive propagator 

To set up the formalism we consider the scalar bubble integral with one massive propaga- 
tor, as shown in Fig. 12 (bubbles 4 and 5 are of this type). The evaluation of the coefficient 
of this bubble integral starts from the identity, 

J dH S^+\l^)S^+\{l - Pf - m") = tdt J {I dl) [I dl] 5{P'' - - t{l\P\I\) , (77) 

so that the formalism follows through as in the massless case, except that t is now frozen at 
the value. 



23 



FIG. 12: Diagram for the calculation of the coefficient of the scalar integral BQ{P;0,m). The 
massive quark is represented by the heavy (red) line. The momenta on the external lines are all 
outgoing. 



The discontinuity of the scalar integral is given by, 



(Idl) [Idl] fP^-m^^ 



mi] V (^1^1^] 
"^^J {i\pm\p\i]- ^^'^ 

Using the standard formula [11, 20] we can express the integrand as a total derivative, and 
perform one of the spinor integrations. 

We obtain, 



where t] is an arbitrary massless momentum. We perform the final integral over (/ dl) by 
inspecting the residues of the pole in the integrand when (/| = [?7|P to find, 

AUP,m) = -iP^-m^)^^^ (82) 

— P^j • (83) 

So, in contrast to the massless case, when applying the spinor integration approach an 
additional rescaling factor of P^/ (P^ — m^) must be applied in order to obtain the coefficient 
of the scalar bubble integral. 
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n 

4 3 

FIG. 13: The box diagram contributing to the rational part in leading colour. 
E. Calculation of the rational terms 

The purely rational terms are calculated using a traditional Feynman diagram approach, 
with tensor integrals handled using Passarino-Veltman reduction [25]. This might seem like 
a retrograde step, ineluctably leading to the algebraic complexity that we have been trying 
to avoid in this amplitude calculation. However two details of our particular case make it 
quite simple. 

First, in the calculation of the rational part we only need to consider diagrams that 
violate the cut-constructibility condition [9]. This states that if n-point integrals {n > 2) 
have at most n — 2 powers of the loop momentum in the numerator of the integrand, and 
the two-point integrals have at most one power of the loop momentum, they will be cut- 
constructible. For our particular calculation the pentagon diagrams are all cut-constructible, 
and, in leading colour, there is only one box diagram, shown in Fig. 13, which is not cut- 
constructible and hence gives a rational part. The calculation of the rational parts from 
lower point diagrams, n = 3,2 does not lead to great algebraic complexity. Second, the 
rational terms are generated by terms in the Passarino-Veltman decomposition that involve 
the metric tensor, g'^'^. Therefore any terms that are not proportional to g'^" may be discarded 
at an early stage of the computation. 

The results for the rational terms turn out to be quite simple. For the { — , +} case we 
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find, 



^^'^(^q ' ^^Q' ' 



+ 



(16) 



2s 



56 



(1 3> [1 5] [1 2] 



(13) [15] / 5234 [12] ^ 
.1(2 + 3 + 4)11] W4|f2 + 3)11] V 



(1|(2 + 3)|1](4|(2 + 3)|1] .23(11(2 + 3 + 4)11] V(4|(2 + 3)|l] 

flip } , (84) 



[24] /[25] (34) [45] N 
S234 V[2 3] (4|(2 + 3)|4]/ 



where the flip operation is deflned in Eq. (27). The contribution for the opposite spin labels 
on the heavy quark line is simply related to this one, 



V' 2q, 3J, 4+, 5+, 67) = R'^il^, 3J, 2^, 4+, 5+, 67; 
For the {— , — } spin labels the result is, 

^''(1."'2o>3q>4+ 5?,67) = -^(-i)4-''(l-,2^,3Q,4+,5+,6,- 



(85) 



+ Af- 



[16) 



2s 



56 



[15] 



(l|2 + 3 + 4|l]s23 



_ ((1|2|1]-(1|3|1]) / [15] [45]x 
(1|2 + 3|1] V(4|2 + 3|l] S23 ^ 



(4|2 + 3|1] 
flip 



where the flip symmetry applies to all terms inside the curly brackets {. . .} but not the 

pref actor J\f , defined in Eq. (14). The final combination of spin labels is obtained by 

symmetry, 

^''(1."' 2J, 3J, 4+ 5+ 67) = -flip (i?^=(l-, 23, 3q, 4+ 5+, 67)) , (87) 
which we note is equivalent to replacing J\f by A/++ in Eq. (86). 



V. THE RESULTS FOR PRIMITIVE AMPLITUDE Af 

The subleading colour primitive amplitude is shown in Fig. 4. Fortunately with our choice 
of massive spinors, Eq. (10), the entire result for all spin labels of the massive quark line 
can be obtained from the one-loop results of Bern, Dixon and Kosower in Eqs. (12.10,12.11) 
of ref. [18], after some replacements and manipulation. The calculation reported in ref. [18] 
was for zero mass quarks and leptons with the following helicity assignments, 

^6 (Ij, 2^, 3+47,57,6+). (88) 
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According to Fig. 2 our standard labelling of the graphs (for one of the amphtudes whose 
spin labelling corresponds to a non-zero amplitude in the massless case) is, 



(89) 



To establish the correspondence between our massive amplitudes and the massless ones 
of ref. [18], we first consider the tree graphs. In our notation the tree amplitudes are given 
by Eqs. (26,28). Note that we must perform the interchanges (1 3), () ■<->■ []), to compare 
with Eq. (12.9) of ref. [18]. As a consequence of our Eq. (12) the massive {— ,+} tree- 
level amplitude, Eq. (26), calculated from diagrams of Fig. 2 and expressed in terms of the 
massless vectors, ki, is identical to the massless result presented in Eq. (12.9) of ref. [18] 
after performing the above interchange. 

This same transformation, flipping the sign of all helicities and interchanging 1 and 3, 
can be used at one-loop level to obtain the bulk of the results for the massive theory. This 
is true for the upper two diagrams of Fig. 4, because the massive quark enters only in the 
form of the heavy quark current see, Eqs. (12,13). The primitive amplitude shown in Fig. 4 
can be split into divergent {V^^) and finite {F^^) pieces as follows 



^6 



(90) 



After performing the interchanges to reduce it to our notation, the singular part of the 
one- loop amplitude as reported by Bern, Dixon and Kosower in ref. [18], Eq. (12.10) is. 



V'\l, 2, 3, 4) = T/''°"(l, 2, 3, 4) + yvertex^^^ 2, 3, 4) 



(91) 



where the two contributions correspond to the upper and lower row of Fig. 4 in the massless 
theory. 



y''°^(l,2,3,4) 

^vertex^^^2,3,4) 



hi 












7i 





3 

3 
2t 



/i 



Sl4 
2 



-4 



'■S23 



7 

2 ' 



(92) 
(93) 



where s, 



2ki ■ kj. 1/™''*"^ = 



in Eq. (93) is the complete correction to external vertex for a 



massless line in the FDH scheme. For the massive case this must be replaced by the vertex 
correction for a massive line, i.e. the result for the lower two diagrams of Fig. 4. The result 
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IS, 

yvertex ^ yvertex ^ yvertex + l^ln^ ^ (95) 

where (3 and /3± are given in Eq. (7), x = — and, 

G'(x) = -2 Li2(-x) -21nxln(l + x) + -ln2(a;) , Li2(x) = - / — ln(l - . (96) 

2 6 Jo 2; 

The self energy corrections on the external massive lines will be accounted for separately in 
association with the wave function renormalization. This concludes our description of the 
divergent parts and the lower two graphs of Fig. 4. 

We now turn to the finite parts of the massive primitive amplitudes shown in the upper 
part of Fig. 4. Because of Eq. (12), the results for the {+, — } and {— , +} amplitudes from 
the upper row of Fig. 4 in the massive theory, expressed in terms of the lightlike momenta 
ki, are given by the results in the massless theory. 

The only remaining issue is whether we can also obtain the massive {— , — } and {+, +} 
amplitudes from the massless results. Note that in ref. [18] the finite parts of the massless 
amplitudes are written in terms of certain symmetry operations in order to make the am- 
plitudes more compact. As a first step we write out the amplitudes explicitly. With this 
result in hand we want to address the issue of whether the results obtained with an external 
fermionic current, Eq. (12) (i.e. the massless one-loop amplitude), can be used to obtain the 
results with the fermionic currents of Eq. (13). Thus expressed in our notation, the massless 
amplitude for helicity choice (2j, 3~) must contain one (3| and one |2]. All other dependence 
on k2 or can only enter in the combination k2 + k^ which can be eliminated by momentum 
conservation. After the amplitude has been recast in this form, we can obtain the required 
result by replacing the current of Eq. (12) with the current of Eq. (13). 

An example may help to clarify the procedure. We shall consider a particular box which 
contributes to the one-loop amplitude, 

^^'(17, 2^^ 3^^ 4+, 5+, 67) = c^4|i|23(^3, ^2) ^0(^4, A;i,P23; 0, 0, 0, 0) + . . . (97) 

The result for this box coefficient in the massless helicity-conserving amplitude, (adapted 
from the first line of Eq. (12.11) of ref. [18]) is, 

, . , ( (13)^[4 5]^ [l|(2 + 3)|6)^[2|(l + 3)|4)n 

'^'^'"^"'^^ " '^^^(32)[56][4|(2 + 3)|1) ' [3 2] (5 6) [1|(2 + 3)|4)^ ) ' ^''^ 
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p 



V 



FIG. 14: Fermion loop contribution to the gluon vacuum polarization. 
This can be rewritten in a form which makes the |3), |2] and k2 + structure manifest, 



C^4|l|23(~) 



-Sl4 



;i 3) (1|(2 + 3)|2][45]' [1|(2 + 3)|6)'(3|(2 + 3)(1 + 2 + 3)|4)[2|(1 + 2 + 3)|4) 



S23 [5 6][4|(2 + 3)|1) 



S23 (5 6) [1|(2 + 3)|4)^ 



(99) 



Replacing the { — , +}-current, Eq. (16), with the { — , — }-current, Eq. (17), we obtain the 
result for the coefficient of this box in the massive amplitude labelled { — , — } 



'^4|l|23(-j 



X 



;i2)(13)[4 5f 



(3 2) [56][4|(2 + 3)|1) 



[1|(2 + 3)|6)^[2|(1 + 3)|4)[3|(1 + 2)|4 ) 
[3 2] (5 6)[1|(2 + 3)|4)^ 



(100) 



Carrying out this operation for all the terms in Eq. (12.11) of ref. [18] we obtain the complete 
massive amplitude for the {— , — } and {+, +} spin labellings. 



VI. RESULTS FOR THE PRIMITIVE AMPLITUDES A^/ AND 



hf 



The unrenormalized contribution of the fermion loop diagram, shown in Fig. 14, for a 
quark of mass m is. 



where (not including the minus sign for a fermion loop). 



n(/,m^) = 2T^l(4 



^ , 4x(l - x) 
ax- 







1 -x(l -a;)^ 

and Tr = |. For nif quarks which can be considered massless this becomes, 

-2 1 / y 10- 



3 e V —p^ 



9 



+ 0(e). 



(101) 



(102) 



(103) 
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Thus in our notation, (c.f. Eq. (31)) the result for the unrenormalized fermion-loop primitive 
is, 

'21 / /i^ Y 10" 



^6^(1^> 2^, 3^, 4q,; 5^ , 6^ ) — AQ^{lg, 2q, 3^, 4^; 5^ , 6^ ) 
where A^^^^ are given in Eqs.(26, 28). 



3 e \-S23 y ~^ 9 



(104) 



VII. RENORMALIZATION 

The amplitudes presented so far are bare amplitudes, which require ultraviolet renormal- 
ization. The renormalization scheme is sUghtly more complicated in the presence of massive 
particles [29] so we specify it in detail here. The requirements for our renormalization scheme 
are, 

• The decoupling of heavy quarks should be manifest. 

• The evolution equations for the running coupling and for the parton distribution func- 
tions should be the same as the equations in the theory without the heavy quark. 
Both the strong coupling and the parton distribution functions should run with the 
coefficients appropriate for the MS scheme in the absence of the massive particles. 

• The mass parameter should correspond to a pole mass. 

These three requirements completely specify the renormalization scheme. If the diagram 
in question contains no heavy internal loops of heavy particles we use the MS scheme. If 
on the other hand the diagram contains heavy loops we will perform subtraction at zero 
momentum, p = 0. The resultant renormalized Green's function will be a function of /w? 
and hence exhibit decouphng as the mass of the heavy quarks becomes large. 

The full renormalized amplitude, A^.^^ is obtained by adding an overall counterterm, 

N,A^.^ = N,A^,^ 

(105) 



We will now justify the contributions in the counterterm term by term. The first term in 
Eq. (105) is the normal MS renormahzation of the coupling constant, which includes the 
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renormalization of the nu loops of massless fermions. The second term in Eq. (105) deals 
with the case where we have rihf heavy fermions. Renormalizing Eq. (102) at = to 
ensure decoupling of the heavy fermions as m becomes large, we obtain for each of the ^hf 
heavy fermions, 

^^(p^m^) = -4 / dx x{l - x) In (l - x{l - x)^) . (106) 
Jo ^ m / 

In the high mass limit simplifies to, 

U^{s,m')^^—^+0{—^). (107) 
15 

Thus after renormalization the contribution of the heavy quark is given by (c.f. Eq. (31)) 

4^(l+,2j,35,4,;5+,67) = 71^(1^ 2j, 3j, 4,; 5+ 67)0^(523, m^) . (108) 

We must also perform a finite renormalization of the gauge coupling [27] to translate from 
the FDH coupling to the normal MS coupling, 

«™" = «r(^l + f . (109) 

This explains the third term in Eq. (105). The last term in Eq. (105) represents the wave 
function renormalization for the two external massive fermions, calculated in Appendix D. 
In the FDH scheme we have from Eq. (DIO), 



1 - g'^crCp 



- + 31n('4l+5 



+ 0{g\e), (110) 



independent of the gauge-fixing parameter in any covariant gauge. 

Since our calculation is performed in the four-dimensional helicity scheme there is a 
further finite renormalization [27] required to arrive at the 't Hooft-Veltman scheme. We 
shall work consistently in the FDH scheme, so this will not be required. 



VIII. IMPLEMENTATION INTO MCFM 



The one-loop matrix elements, computed using the methods described above, have been 
included in a full next-to-leading order calculation of the WQQ process. The amphtude for 
the lowest order process is given in Eqs. (26,28). In order to complete the NLO calculation 
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FIG. 15: Feynman graphs appearing in the calculation of the real radiation amplitudes. 

the Born level amplitude and the one-loop amplitude must be supplemented with results 
for the real radiation diagrams and a method for cancelling infrared singularities between 
the two contributions. The tree-level real radiation process has been computed using the 
diagrams shown in Fig. 15, adopting the same choice of massive spinors as used in the 
virtual contribution. Infrared singularities are handled using the subtraction method [32] 
implemented using the dipole formulation [33] and extended to the case of massive emitters 
and spectators [34]. The full calculation will be made available as part of the the MCFM 
code [30, 31]. 

To provide a point of comparison we have timed the evaluation of the complete one-loop 
contribution, i.e. the interference in Eq. (25), summed over all four possible spin labels on 
the heavy quark line. Using a standard 2.66 GHz machine and compiling the code with 
the -02 flag of gfortran this operation takes 4ms for a single phase space point, which 
is about a factor of sixty slower than evaluating the corresponding massless amplitudes as 
implemented in MCFM. 
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mw = 80.44 GeV 



nib = 4.62 GeV 



mt = 172.6 GeV 



Gf = 1.16639 • 10-5 GeV 



—2 



sin^ e,aj = 0.223 
a = ^Gprn^ sin' 



w 



Vud = Vcs = 0.974 GeV 



Vus = Vcd = 0.227 



PDF set: CTEQ6L1 



LO 



PDF set: CTEQ6M 



NLO 



as{mw + ^rub) = 0.130345 



as{mw + "^rnb) = 0.118298 



TABLE IV: Inputs used for the 



results presented in section VIII. 



A. Checks on the calculation 

A number of checks have been performed at various stages of the calculation. 

1. At the level of the one-loop amplitude, we have cross-checked our calculation with 
a numerical implementation of generalized (/-dimensional unitarity [35]. This check 
confirms the coefficients of all the scalar integrals and rational terms, as well as the 
complete (unrenormalized) amplitude, at any given phase space point. The values of 
the unrenormalized one-loop amplitudes at a specific phase space point are given in 
Appendix E. These numbers may be useful in performing a check of our calculation. 

2. The cancellation of infrared singularities is performed using a slight extension of the 
original dipole formulation in which the extent of each subtraction region is controlled 
by an additional parameter [36, 37]. This parameter also appears in the integrated 
form of the dipole counterterms in such a way that the sum of real and virtual radiation 
does not depend upon its value. We have checked that this independence is indeed 
manifest in our calculation. 

3. We have checked that our final results for the Wbb integrated cross section agree with 
the values reported in the earlier calculations of Refs. [2-4] . 

B. Phenomenology 

For now we present only a limited set of results, focussing mainly on the comparison with 
previous calculations and presenting a simple distribution involving the decay products of 
the W boson. We leave a detailed phenomenological study for a future work. 
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7 TeV 


8 TeV 


14 TeV 


LO 

W+bb 

NLO 


4.456(2) 
8.655(9) 


5.157(2) 
10.58(2) 


9.041(3) 
23.51(3) 


LO 

W'bb 

NLO 


2.588(2) 
5.053(5) 


3.109(1) 
6.353(6) 


6.256(2) 
15.55(2) 



TABLE V: LO and NLO cross sections (in picobarns) for Wbb production at various energies of 
the LHC. Integration errors are shown in parentheses. 

For convenience, we choose the same set of input parameters here as reported in Ref . [3] , 
which are summarized in Table IV. The final state is defined by the following cuts on the 
6-jets, 

> 25 GeV , Iry^l < 2.5 , (111) 

where the jets are identified using the kx clustering algorithm with pseudo-cone size R — 0.7. 
The results presented here are inclusive of the additional jet that may be present at NLO. 

As already discussed previously, although we treat the 6-quark as a massive particle when 
it appears in the final state, we use nif = 5 light flavours in the running of ctg and the PDF 
evolution. This is primarily for comparison with previous work [3]. We note that because of 
the smallness of Vcb and Vub the 6-quark distributions in the initial state make a negligible 
contribution. However because S23 > 4m^ it is more appropriate to have a strong coupling 
constant running with riif = 5 active flavours. 

We first present the cross sections for this process at the LHC, for center-of-mass energies 
a/s — 7, S and 14 TeV. The results at LO and NLO are shown in Table V where we have 
used a scale choice fJ.R = fJ^F = mw + Sm^ throughout, again to facilitate comparison with 
Ref. [3]. For these results no cuts arc applied on the decay products of the W boson and 
the corresponding branching ratio is removed, so that the reported cross sections are for a 
W that does not decay. In this way, one sees from the final column {^/s — 14 TeV) that 
our results agree with those reported in Ref. [3] at the level of 0.5%. As noted in earlier 
work [38] , at LHC energies the NLO corrections are substantial because of the influence of 
the quark-gluon initial state. 

To illustrate a new quantity that may now be computed at NLO accuracy, in Fig. 16 we 
show the LO and NLO predictions for the quantity i?^™™ which is defined as the separation 
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1 2 3 



min 



FIG. 16: Distribution of the minimum separation between the electron and any jet, Rij, defined in 
Eq. (112) of the text, for W~bb production at the 14 TeV LHC. The distributions are normalized 
to unit area. 

between the charged lepton and the closest jet, 

= min J (r^iepton - Vict)^ + (01epton " 0jet)^ • (112) 
{jets} V 

In this equation the azimuthal angles and pseudorapidities (in the lab-frame) of the lepton 
and jets are denoted by and r] respectively. We see that the effect of the NLO corrections 
on the shape of this distribution is relatively minor, visible only for R^^"^ > 1.5. 

IX. CONCLUSIONS 

We have presented the first computation of the Wbb cross section with massive 6-quarks, 
including the lepton correlations present in the decay of the IV-boson. This calculation 
required knowledge of the one-loop virtual corrections to the qQQqIi process retaining the 
mass for the heavy quarks Q. The one-loop amplitude was obtained using the spinor helicity 
formalism. This method has been extensively used for one-loop calculations with massless 
quarks, but rarely with massive quarks. The calculation required a number of modifications 
of standard techniques to cope with the presence of the mass. 
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Although our results are analytic we have not yet simplified them sufficiently to publish 
them in a journal article. Our results for the one-loop amplitudes will be included in the 
released version of MCFM, which is an appropriate method of publishing such results. Our 

analytic results did lead to a code which is fast and numerically stable. Using this code we 
intend to study the detailed phenomenology of this process in a future publication. 
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Appendix A: Notation for spinor products 

The spinor notation is almost standard, but because of the presence of massive particles, 
it is important to make it explicit. We adopt the following notation for massless spinors, 

(i| =: (i - I = u-{ki), [i\ = u+{ki) . (Al) 

Further the spinor products are defined as, 

(ij) = (i - \ = u-{ki)u+{kj) , 

[ij] ^ {i + = u4h)u4kj) , (A2) 
with ki, kj massless particles. With our convention, 

{ij) [ji] = • kj = Sij 

We also define the spinor strings. 

{t\j\k] ^ {k,-\l^,\kk-) , 
{i\j + k\l] = {k, - \{^j + J^,)\ki-) , 
{i\jk\l) = {ki - \f^jh\ki+) ■ 
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(A3) 



(A4) 



Pa -Pz 




(a) (b) (c) 

FIG. 17: Diagrams for two gluons coupling to a massive quark line. 

For the case of a massless momentum kj we may write, 

mk] = (ij)[jk] , (A5) 

but for the case of a massive momentum, (in our notation the momenta p2 and ps), this 
separation in no longer possible. As a compact notation we therefore write in this case 
(i|jffc], denoting the massive momentum by a bold-face symbol. 
The Schouten identity, 

[i k] [m n] = [i n] [m k] + [i m] [k n] , (A6) 
may be applied to these compound quantities. Thus we have, 

(2|jf/c][mn] = (i|jfn][mA;] + (z|jjm][/cn] . (A7) 



Appendix B: Tree level results 



1. Results for A{lg, 2q, 3q, 4J 

As an example of the use of the massive spinors employed in the calculation we consider 
the tree-level two-quark two-gluon amplitude, A{lg, 2q, 3q, 4^). The fermions, with momenta 
P2 and have a common mass m. The momenta are all outgoing so that ki+p2+P3 + k4 = 0. 
The amplitude is obtained by summing the contributions of the three Feynman diagrams, 
M("), M(*) and M^^) shown in Fig. 17, 



A 



tree 



-.9 



Normalizing the colour matrices as follows, 

Tr T^T^ = , 



(Bl) 



(B2) 
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an explicit calculation yields the following contributions, 



(B3) 



^(5) ^ £i:4i^(p^)^^^(p^) 



ki ■ ki ■ k4^ 



^'^U{p3)^iv{p2) + ^^-p-u{p3)^4v{p2) , (B4) 



Pl3 ~ 



(B5) 



where pia = /ci +P3 and P34 = pa + k^. Prom Eq. (10) we write the massive spinors as, 

u{ps) = Nsu{k2M + m) , (B6) 
v{p2) = A^2(y2-m)^;(A;3) , (B7) 

and then reorganize the calculation by using the identities [39] , 

(^3 + m)^4 = (j^u - |!f4 + m)^4 

= (2p34 - k^) ■ 64 - ^[1^4, ^4 - ^4(^34 - m) , (B8) 

(-i/12 + = -{2p2 + ki) ■ ei - ^[^1, /i] + /i(^2 + m) . (B9) 



Inserting these identities we may write new forms for M^"'^ and M*^'^^ 
M« = N2N3u{k2) r(') (^2 - m) t;(A;3) , i = {a, b, c) 



p(a) 

r(c) 



I [2^3-64 -^[^4,^4] 

{[2p3-ei-^[^i,?i] 



2p3 • k4 

1 



- 2p2 • ei - ^[^1, /i] - /4/1 } , 

- 2p2 • £4 - ^[^4, ^4] - ^Wa^ ■ 



(BIO) 
(Bll) 

(B12) 



2p3 • h 

These formulae have the advantage that the m dependence is corralled in a single place, 
in Eq. (BIO). Moreover the overall form is simple and the last term in Eqs. (Bll, B12) 
does not contain the massive propagator. F^^^ can be read off from Eq. (B4). Inserting the 
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appropriate polarization vectors we find, 



-iA{lg ,2^,3q 
-iA{lg ,2q,3+ 

-iA{lg , 2t, 3q 
,2^,3+ 

2. Results for A(l5,2g,3g) 



[2 4f (1|;^|4] 
[14][2 3](4|^|4] ' 

(12)^(1|^|4] 
(14)(2 3)(4|^|4] ' 

m (13)^[34]^ 



4+) = 2(/3+-/3_) 
4+) = 2(/3+ - /3_ 

4;) = , 

4,-) = 0, 

(2 3)(4|^|4] 
^ m/j- (14)^ 
[2 3](4|^|4]- 



[2 3] (14)[14](4|^|4] ' 

m (13) '[3 1]' 
'(2 3) (14)[14](4|^|4] ' 



(B13) 



The tree-level results for the simple qqg amplitudes, stripped of overall colour and coupling 
constant factors, are well known: 



-a(i-,2+,3;) 

-a(l- 2+ 3t) 



(13)^ 
(12) 

[2 3]' 



(B14) 
(B15) 



3. Results for A{lq, 2,, 3^, 4;-, 5/) 



The tree-level amplitudes for qqWg are also rather simple. Removing the colour and 
coupling constants as normal we have. 



-z^(l-,2+,3;,4±,5p 
-z^(l-,2+,3+,5±,4-) 



[2 4]^ 



[13] [2 3] [4 5] ' 
(13) (2 3) (4 5) • 



(B16) 
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Appendix C: Scalar integrals 



We work in the Bjorken-Drell metric so that I — Iq — If — I2 — l^- The definition of the 
integrals is as follows, 



1 



BQ{pi]mi,m2) = — — / d^l jT^ 2\((i , ^2 2\ 

m^rvJ - ml){{l ^ pxf - ml) 

A-D 



Co{pi,p2; mi, 7712, 7713) 

X J d^l 



in 2 rr 



(CI) 



{P - ml){{l + piY - ml){{l + pi +P2y - ml) ' 



Do{pi,p2,P3;mi,m2,m3,m4) 
X / d^/ 



A-D 



D_ 

m 2 rr 



- ml){{l+pi)^ - mi)((Z +pi +p2)2 - +pi + p2 - ml) ' 

We have removed the overall constant which occurs in £)-dimensional integrals 
r2(l-e)r(l + e) 



rr 



r(l - 2e) r(i - e) 



1 + 



+ 0{e') . (C2) 



Appendix D: Fermionic self energy 



Introducing the renormalization of the bare parameters, mo — Z^m and Qo = v ZqQ 
we may write the renormalized inverse propagator for a heavy quark of momentum p as. 



iV^{p,m;g) = Zq - m -T.{p,m;g) -m{Z^- 1) +0{g^) . 



(Dl) 



—iL{p,m;g) is the contribution of the one-loop heavy quark self-energy graph, which prior 
to renormalization has the form, 

E(p,m;^) = -/Ci.cr[x(p2) - m) + mY {p')^^ . (D2) 

By direct calculation of the Feynman diagram in an arbitrary covariant gauge specified by 
gauge fixing parameter A, 



X{p') 



Y{p') 



2(1 - 6e){Boip;0,m) + 5i(p;0,m)) - (1 - X){Bo{p;0,m) + {p' - m')B'^{p;^,m)) 



2Bi{p,Q,m){l-6e)-2Bo{p;0,m) - (1 - A)(p' - m')5^,(p; 0, m)) 



(D3) 
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where 5 — in the FDH scheme, 5 — 1 in the conventional dimensional regularization 
scheme, and Cp — ■ '^^^ integrals Bq and Bi are defined as. 



{So(p;0,m),Si(p;0,m)p''} 



ITT 2 rr 



/2 ((/+p)2_m2) ' 

and B[ is the derivative of the form factor Bi with respect to p^. 

Taking the limit — rv? before the limit e ^ 0, we have the following results. 



(D4) 



5o(p;0,m)|p2=^2 = (-^ 

(11^ 

-Bi(p;0,m)|p2=„2 = — 



2. .1 



- + 2 

e 



Si(p;0,m)|p2=^2 



2e 2J 

2. .1 



2m2 
1 



- + 2 
e 



2m2 



(D5) 



The mass renormalization is fixed by the condition that the inverse propagator vanish on 
shell. 



Z^^\^g'CFY[m^)^\-cvg'Cp - + 31n(^)+5-5 ^0{g\e) 
After mass renormalization the result for the inverse propagator becomes, 

-ir^(p, m; g) = Zq [{^ - m)(l + /C^crX(p2)) + m C^cr(y(p') - ^(m^)) 
Renormalizing the wave function at the point p^ — we find. 



1 - g^CpCT 



dp'^ 



(D6) 



(D7) 



(D8) 



By explicit calculation we have that. 



1 , / /i' 



+ In 



2m 



,dY{p^ 



dp^ 



+ (3 - 5) - (1 - A) - + In 



1 , / /i' 



+ 2 +0(e) 



+ 2 In 



+ 2 + (1 - A) I - + In 



/i 



+ 2 +0(e) 



(D9) 



The final result for the wave function renormalization is independent of the gauge fixing 
parameter. A, 



1 - g^crCp 



- + 31n( ^ 1 + 5-5 



+ Oig\e). 



(DIO) 



The agreement between Eq. (D6) and Eq. (DIO) is fortuitous because the former contains 
only ultraviolet poles, whereas the latter contains both ultraviolet and infrared poles. 



41 



Appendix E: Numerical evaluation 



For the convenience of the reader we present numerical results for the amplitudes at a 
particular phase space point. The results we present contain no ultraviolet renormalization 
and the self energy corrections on the two massive external legs have not been included. 
(The self energy corrections on the massless external legs vanish). 

The phase space point is specified by six momenta satisfying overall momentum conser- 
vation ki + p2 + P3 + + + kg — and /cf — 0, P2 — pI — The massive momenta 
P2 and p3 are defined in Eqs. (6,7). For the numerical results the massless momenta ki are 
taken to be, 

ki — ^(— 1, -|-sin^, -|-cos^sin0, -|-cos^cos0), 

k2 = |(1, 1,0,0), 

ka — —(1, cos (7, sin (7, 0), 

A;4 = ^(— 1, — sin^, — cos^sin^, — cos^cos^), 
k^ = ^(1, cos p cos (7, cos p sin o", sin p), 

ke = —ki — k2 - ks — ki - k^ , (El) 

where = 7r/4, = 7r/6, p = 7r/3, cosu = —7/19, m = 1 GeV and p. = 6 GeV. Note that the 
energies of ki and k^ are negative and kf — 0. As usual /i also denotes the scale which is 
used to carry the dimensionality of the £)-dimensional integrals. With m — 1 GeV we have 
/3 ~ 0.38397382 ... so that /3 differs substantially from the massless hmit, /3 — 1. 

The results for the primitive amplitudes are presented in Tables VI and VII, where we have 
divided the 1-loop amplitudes by their tree-level counterparts in order to remove the overall 
ambiguity in the phase [40] . The individual amplitudes labelled by the particular spin labels 
are dependent on our choice of the spinor wave functions, Eq. (10). Results independent 
of this convention may be obtained by summing the squares of the four amplitudes. The 
results in Tables VI and VII have been checked by an independent program. 
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